Superdiffusive transport by multivalent molecular walkers moving under load 
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We introduce a model for translational molecular motors to demonstrate that a multivalent catalytic walker 
with flexible, uncoordinated legs can transform the free energy of surface-bound substrate sites into mechanical 
work and undergo biased, superdiffusive motion, even in opposition to an external load force. The walker in 
the model lacks any inherent orientation of body or track, and its legs have no chemomechanical coupling other 
than the passive constraint imposed by their connection to a common body. Yet, under appropriate kinetic 
conditions the walker's motion is biased in the direction of unvisited sites, which allows the walker to move 
nearly ballistically away from the origin as long as a local supply of unmodified substrate sites is available. 
The multivalent random walker model is mathematically formulated as a continuous-time Markov process and 
is studied numerically. We use Monte Carlo simulations to generate ensemble estimates of the mean squared 
displacement and mean work done for this non-ergodic system. Our results show that a residence time bias 
between visited and unvisited sites leads to superdiffusive motion over significant times and distances. This 
mechanism can be used to adapt any enzyme-substrate system with appropriate kinetics for use as a functional 
chemical implementation of a molecular motor, without the need for structural anisotropy or conformationally 
mediated chemomechanical coordination. 
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I. INTRODUCTION 

Molecular motion at the nanoscale is dominated by ran- 
dom, thermally driven collisions that lead to slow, uncontrol- 
lable diffusive transport. Diffusion of large cargo molecules 
in a crowded cellular environment is so slow that nature 
has evolved sophisticated, specialized molecular machines to 
transport cargo directionally at superdiffusive rates. These 
translational molecular motors, such as kinesin, dynein, and 
myosin Q, walk along oriented tracks, consuming chemical 
energy in the form of ATP and converting it to mechanical en- 
ergy that is used to do work against external load forces 0. 
Molecular motors are essential to a cell's ability to control 
the distribution and movement of materials and information 
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FIG. 1 : (color online) A multivalent random walker (MVRW) has a 
rigid body and k flexible legs of length t. A leg can attach to and de- 
tach from fixed chemical sites that are within length £ from the body. 
The enzymatic action of the leg can irreversibly transform a substrate 
site into a product, changing the subsequent binding kinetics for the 
site. As the legs attach to and detach from sites, the walker moves 
over the surface. 
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within it. 

Synthetic nanoscale systems, like the natural cellular sys- 
tems that have inspired their development, also need mech- 
anisms to maintain non-equilibrium distribution of materi- 
als and information [4]. Useful synthetic behavior has been 
demonstrated by combining natural molecular motors and 
components in novel ways (5). However, each species of nat- 
ural molecular motor is highly specialized to its cellular en- 
vironment, chemical fuel source, and polymeric track (e.g., 
microtubules, actin, or DNA). Synthetic systems with differ- 
ent polymers must either adapt existing natural motors, or use 
newly designed compatible synthetic motors. Kinesin and 
other natural motors can be mutated to change the kinetics of 
their motion O, but not to catalyze arbitrary fuel substrates, 
or move over arbitrary tracks without fundamentally alter- 
ing their functionality or efficiency. Natural molecular mo- 
tors rely on non-local conformational changes to couple the 
binding of fuel with the kinetics of track binding (7), and use 
this conformationally-mediated chemomechanical coupling to 
coordinate their processive hand-over-hand walking gait (8). 
These mechanisms make natural motors efficient but also hard 
to mimic in synthetic systems. 

We show that neither oriented tracks, nor rigid walking 
gaits, nor chemomechanical coupling, nor coordinated confor- 
mational changes, are necessary for a molecular walker to do 
ordered mechanical work as a molecular motor. We consider 
the motion of a multivalent walker with a rigid, inert body 
and several flexible, enzymatic legs (Fig.[T]). The legs attach 
to and enzymatically modify surface-bound chemical sites ar- 
rayed as arbitrary 2D patterns and tracks. The legs are chem- 
ically and conformationally uncoupled, other than the passive 
constraint imposed by the connection to a common body. Yet, 
under appropriate kinetic conditions, the walker can be made 
to move directionally and processively even in opposition to 
a force. By modeling and understanding such simple walker 
systems, we learn which chemical and mechanical properties 
of walker-based motors are sufficient for superdiffusive mo- 



tion, and which properties are not necessary. 

Our multivalent random walker (MVRW) model is a gen- 
eral description of the motion of these walkers. It takes the 
form of a continuous time Markov process that describes the 
stochastic motion of a walker as it moves over and modifies 
surface-bound sites. The model is designed so that the Markov 
process can be efficiently simulated, yet still accurately model 
the effect of external forces on the physical motion and chem- 
ical action of the walker's body and legs. 

Through Monte Carlo simulations we show that walkers 
with appropriate kinetics can move superdiffusively in the 
direction of unvisited sites over significant times and dis- 
tances, and can do so while performing a non-trivial amount of 
work against an external load. This effect can be understood 
from the spontaneous emergence of a substrate concentration 
gradient — a boundary between visited product sites and un- 
visited substrate sites. At this boundary, the non-uniform local 
substrate concentration combined with the chemical kinetics 
of the legs and the constraints that the body places on leg mo- 
tion lead to a directional bias away from previously visited 
sites. As the legs irreversibly modify substrates to products, 
they move this bias-inducing boundary further from the ori- 
gin. Hence, as long as a walker stays proximate to the bound- 
ary it moves ballistically away from the origin. 

The MVRW model is inspired by attempts to model 
the motion of synthetic DNA-based molecular walkers, 
called molecular spiders . The spiders use a deoxyri- 
bozyme/nucleotide ifTOli enzyme/substrate chemical imple- 
mentation. Molecular spiders have been observed to walk 
processively in 3D environments (9|, and move direction- 
ally over 2D nanoscale tracks ifTTIl . Abstract models of ID 
spider motion [Q21 El have been shown to exhibit tran- 
sient superdiffusive behavior under appropriate kinetic con- 
ditions |[T4ll . Other ID models have studied the effect of force 
on dissociation rates using the kinetic rates specific to deoxyri- 
bozymes [ 151 [T6l . 

The MVRW model shows that mechanisms for designing 
molecular motors exist without the need for chemomechanical 
coupling, conformational coordination, rigid walking gaits, or 
inherent orientation. Multivalent random walkers, like natu- 
ral molecular motors, are Brownian ratchets ifTTl that rectify 
random molecular motion into ordered work and directional 
transport. Both systems achieve this rectification by utilizing 
the chemical free energy of a substrate fuel. However, the 
mechanisms by which MVRWs do this are significantly dif- 
ferent from natural motors. Unlike kinesin and other natural 
motors, MVRWs move over arbitrarily arranged 2D tracks, 
and are able to do so without inherent orientation or structural 
asymmetry. The gaits of a molecular walker are uncoordi- 
nated and acyclic, yet the irreversible modification of surface 
sites causes an emergent asymmetry in local substrate con- 
centrations that is able to bias the motion of walkers, allow- 
ing them to move directionally along prescriptive landscapes. 
The structural and chemical simplicity of MVRWs is one of 
their most important properties as it means that the conceptual 
functionality of a molecular spider is independent of the spe- 
cific enzyme/substrate system used in their implementation. 
Hence, multivalent random walkers provide a different per- 



spective for better understanding what structures, properties, 
and mechanisms are minimally necessary to turn a molecular 
walker into a molecular motor. 



II. THE MULTIVALENT RANDOM WALKER MODEL 

We model the motion of a multivalent random walker as a 
continuous-time, discrete- state Markov process, where each 
state transition corresponds to a chemical reaction between a 
leg and a surface-bound site. The model conceptually sep- 
arates the timescales of the relatively slow leg- site interac- 
tions from the much faster physical (mechanical) motion of 
the walker's body and legs. In this way, only chemical re- 
actions correspond to state transitions in the Markov process, 
and the state space remains discrete. The Markov property is 
ensured by assuming that the physical motion of the position 
of the body and unattached legs comes to an equilibrium in 
between successive chemical reactions. 

The MVRW model is 2D and consists of a walker and an 
environment of surface-bound sites. The walker has a rigid, 
point-like body which serves as the attachment point for k 
flexibly tethered legs, each with maximum length L The en- 
vironment is described by a set S C M 2 of immobile chemical 
sites. All sites are initially substrates, but they can be trans- 
formed into products by the enzymatic action of a leg. 

The state of the Markov process needs only to describe the 
state of the reactive chemicals in the system, i.e, the species 
at each chemical site and the chemical state of each leg. The 
state of the surface is defined by the set P of sites that have 
been transformed to products; all other sites in S are consid- 
ered to be substrates. A walker leg is either attached to a site 
in S or is detached. No two legs may be attached to the same 
site. The state of the walker is succinctly represented by the 
number of detached legs < d < k and the set A C S of at- 
tached sites. Thus, any state CO can be described compactly by 
a triple CO = (P,d,A), and we let £1 be the set of all potential 
states. 



A. Chemical kinetics and state transitions 

The kinetics of a leg (L) binding to substrate (S) and prod- 
uct (P) sites is controlled by five reaction rates ll43l 

L + S^LS ^>L + P^LP. (1) 

k s k p 

The reactions of Eq. [T] correspond to three types of func- 
tional motion for walker legs, association (binding), dissoci- 
ation (unbinding), and catalysis. Each of these actions corre- 
sponds to a transition in the Markov process. In the MVRW 
model, we assume that rates for the unimolecular dissociation 
and catalysis reactions for a leg-site complex are not chemo- 
mechanically coupled to the conformations or positions of the 
body and other legs. Thus, as in the unimolecular stochas- 
tic kinetic models of Gillespie lfT8l . each individual LS or LP 
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FIG. 2: (color online) The feasible body positions, J 7 , as determined 
by the attached leg constraints are shown in yellow. Any site at most 
distance £ from J 7 is a feasible site (green region). 



complex will dissociate or undergo catalysis according to the 
rates fcT , kp , and k cat . 

The bimolecular association reactions are more compli- 
cated as their propensity depends on the likelihood of the leg 
being proximate to the chemical site in order that it may bind. 
This likelihood, in turn, depends on the position of the body 
and the unattached legs. 



B. The body's equilibrium position 

In between reaction events, the walker's body is subject 
to thermally driven constrained diffusion, and we assume it 
comes to a physical equilibrium distribution B. This assump- 
tion ensures that the Markov property holds for the MVRW 
stochastic process because the high-frequency physical mo- 
tion of the walker quickly removes any conformational mem- 
ory of previous chemical states, and the body distribution only 
depends on the current state CO = (P,d,A). 

We assume that the only coupling between the attached legs 
and the body is that the body is constrained to stay within 
distance £ from each attached site. Hence, P [B = p] = for 
position p if there is any attached leg site s G A such that ||p — 
s|| > £. We call all values of p that satisfy ||p — s|| < £ for all 
s G A the feasible body positions, T, illustrated in Fig. [2] 

At equilibrium, B is a Boltzmann distribution over the fea- 
sible positions p G T according to the energy t/(p) at each 
position, 



P[B = p]= jPB (p) 



-pu(j>) 



J>-/Mp)dp' 



(2) 



In Eq.|2j j3 = 1/ksT, where kg is Boltzmann's constant and T 
is absolute temperature. The energy U (p) necessarily depends 
on the mechanical properties of the walker legs and their teth- 
ers. To study the case of minimal chemomechanical coupling, 



we model the energy as uniform over any feasible position, 
but infinite for infeasible positions, 







p G T 

otherwise 



(3) 



C. Leg-site binding kinetics 

The bimolecular kinetics of leg-site binding is controlled 
by two factors, (I) a second-order process by which the leg and 
site come into contact, and (II) a first-order process wherein 
the leg and site undergo conformational changes to move to 
a strongly bound state 1 19 ]. We consider the case when the 
legs are short enough, and the conformational changes lead- 
ing to binding are slow enough that factor II is limiting. In this 
case, an unattached leg undergoing constrained diffusion has 
the opportunity to interact many times with the local feasible 
sites before it finally binds strongly enough to be considered 
attached. The overall rate of a leg reacting with any feasi- 
ble site is then proportional to the number of feasible sites in 
its proximity. Equivalently, from the perspective of a feasible 
site, the probability that it reacts with the leg is independent 
of the number of other feasible sites in the local environment. 
For any site s and body position b we define a feasibility func- 
tion, 
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1 ||s-b||<£ 
otherwise 



(4) 



Then from position b, an unattached leg binds to site s, with 
species 7r(s) G {S,P}, with rate 



n.M=# w /.(b). 



(5) 



Now, we take into account that the body is not at a single 
position b, but in an equilibrium distribution B over positions, 
and we integrate Eq.|5]to obtain 

m(s) = k+ {s) J^p B (*>)fs(y>)db. (6) 

Any site with non-zero rate of attachment is called a feasible 
site', the region of feasible sites is shown in Fig. [2] 



D. The effect of force on walkers 

The MVRW model can also capture the effect of forces on 
the walker body. This is an advantage of modeling the body's 
position as a Boltzmann distribution determined by the energy 
of the walker at each feasible position. Under the effect of a 
conservative load force f, the energy of position p is 



t/f(p) = £/o(p)-f-(p-Po), 



(7) 



where po can be any reference point. The adjusted energy 
Uf gives a new equilibrium distribution with probability mass 
shifted in the direction of the applied force. The effect of 
force on the body's equilibrium position and on the attach- 
ment propensity for each of the feasible sites is illustrated in 
Fig.0 



/ = O.OpN 



/ = 0.5pN 




FIG. 3: (color online) The equilibrium body position probability den- 
sity for a walker under several different load forces. Warmer colors 
represent increasing probability. The body is drawn at its mean equi- 
librium position, (B). The region of feasible sites is illustrated in 
green, but not all sites are equally likely for attachment. The color 
and size of feasible sites indicate the effective attachment rate as de- 
termined by the body position distribution B. 



III. SIMULATION AND ANALYSIS METHODS 

The MVRW model is a continuous-time Markov process 
(CTMP) with discrete states CO = (P,d,A) G £1 Given all 
relevant parameters, the walker CTMP defines certain ran- 
dom variables {X(t)} t >o over the state space £1. This sin- 
gle Markov process simultaneously describes the fast phys- 
ical motion of the walker body and legs under an external 
load, as well as the slower discrete chemical state changes 
of the legs binding and modifying the surface sites. The 
MVRW simulation process uses a kinetic Monte Carlo al- 
gorithm (20j to sample from the non- equilibrium behav- 
ior of the walker's chemical motion, but each chemical step 
in this MVRW Markov Chain requires using the Metropolis- 
Hastings algorithm to sample the equilibrium behavior of the 
body's physical motion. Hence, by separating the timescales 
of the chemical from the physical behavior of the system, we 
can take advantage of both equilibrium and non-equilibrium 
Markov chain Monte Carlo techniques, using each technique 
where it is most applicable to the dynamics of the walkers. 

Our kinetic Monte Carlo (KMC) algorithm for the MVRW 
process produces samples x(t) for t G [0,£ max ], such that at 
each time, x(t) is a sample of a random variable X(t). With an 
ensemble of TV samples of the Markov process, we measure 
and report various properties of the system state at linearly 
and logarithmically spaced time points t G [0, t max ] . The sim- 
ulation algorithm is described in detail in previous work l22l . 



In order to compute the transition rates for the association 
reactions, we use the Metropolis-Hastings algorithm I23H251 
to sample from the body's equilibrium distribution B and use 
these samples for Monte Carlo integration of Eq. [6] Impor- 
tantly, the Metropolis-Hastings algorithm is able to sample 
from B using only the energy from Eq. [7] without having to 
compute the partition function of the Boltzmann distribution 
ofEq.g 

There are computational advantages to modeling the dis- 
tribution B at equilibrium. When the forces on the walker 
are conservative, B is translationally invariant, and depends 
only on the relative locations of the attached legs, and not 
on the whole system state CO = (P,d,A). When the walk- 
ers move over regular lattices, there are only a finite num- 
ber of potential leg attachment gaits, and their corresponding 
attachment propensities can be precomputed, eliminating the 
need to run Metropolis-Hastings at every KMC iteration. This 
makes KMC simulation tractable for long times and large val- 
ues of N. 



A. Random number generation 

In all Monte Carlo methods the fundamental source of 
stochasticity derives from a deterministic pseudo-random 
number generator. The statistical properties of the pseudo- 
random number source are critically important to the correct- 
ness of model predictions [26]. To allow parallel computa- 
tion and preserve mathematical guarantees of random number 
generator quality, we use the leapfrogging method to gener- 
ate N parallel random number streams from a single master 
stream l27l . Hence, only a single random seed is needed to 
compute all N KMC traces for each set of model parameters 
studied. 



B. Mean squared displacement and diffusion 

In single-particle tracking, the stochastic motion of indi- 
vidual molecules is frequently analyzed in terms of the mean 
squared displacement (MSD) (281 . The MSD is the vari- 
ance in the displacement, Var (||p(f)||) = (||p(7)|| 2 ). For any 
diffusive process (i.e., an unbiased random walk) the MSD 
will scale linearly with time. Anomalous diffusion l29l is 
characterized by the MSD scaling as some non-linear power 
0<a<2, 



IIpWII 2 }-^ 



a = stationary 

< a < 1 subdiffusive 
a = 1 diffusive 

1 < a < 2 superdiffusive 

a = 2 ballistic or linear 



(8) 



MSD can either be computed as a temporal average (over dif- 
ferent St values for a single walker trajectory) or an ensemble 
average (over absolute t for an ensemble of trajectories from 
identical walker systems). Many biological systems are (or 
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FIG. 4: (color online) A snapshot several hundred steps into a 
MVRW simulation. The surface track for walkers in this set of sim- 
ulations consists of a semi-infinite strip of substrate sites 3-wide. 
Shown are the circular constraints imposed by the attached legs, and 
the probability density for the body's equilibrium position (as a 
heat map). The walker has one unattached leg, and the relative rate 
at which it would attach to each site is shown by the size and color 
of the site. 



are at least assumed to be) ergodic in the sense that the mo- 
tion of a walker is independent of its absolute position on the 
track and does not depend on its previous motion over a re- 
gion of that track. Under the assumption of ergodicity the 
temporal and ensemble MSD are equivalent (assuming suffi- 
cient measurement resolution), but when a non-ergodic sys- 
tem is analyzed, only the ensemble average is meaningful for 
use in characterizing anomalous diffusion J30J [31) . MVRWs 
are a nonergodic system because they irreversibly modify the 
track as they move over it. Thus, the motion of the walker 
depends on its absolute position on the track and specifically 
on whether the local sites are products or substrates. Hence, 
only the ensemble MSD can be used to study MVRWs. 



IV. RESULTS 

By itself a multivalent random walker is just a rather unso- 
phisticated multivalent enzyme, but when paired with an ap- 
propriately designed nanoscale track of substrates it becomes 
a molecular transport device, able to move superdiffusively 
even under the influence of an external load force. Using 
KMC simulations we studied the motion of MVRWs mov- 
ing over a 3-wide semi-infinite track of substrates (Fig. [4]). 
The relevant simulation parameters are summarized in Ta- 
ble [I] The walker starts with a single leg attached to the mid- 
dle leftmost site. The remaining legs quickly attach, and the 
walker begins to move over the surface. From this initial po- 
sition, the lack of substrates to the left means the walker can 
only move in the -fx direction, so we apply a force in the 
— x direction to oppose the walker's motion. If the force ap- 
plied to the walker is f = (f x ,fy), we let f y = 0, and write 
/ = —f x as a scalar for the magnitude of the force in the — x 
direction. We limit / < 4.0pN as larger forces result in in- 
significant motion under the parameters of Table [I] The up- 
per bound of / = 4.0 pN is near the maximum force a DNA- 
based realization of a MVRW could a priori be expected to 
move against, as the stall force for kinesin is approximately 
5 — 8pN (32), and the dissociation force for double- stranded 
DNA is < 12pN (33). 

Under the kinetics from Table [TJ & cat serves a special role as 



TABLE I: Model parameters used for simulations. 



Parameter Description 


Symbol 


Value 


Number of legs 


k 


4 


Leg length 


£ 


12.5 nm 


Track width 


- 


3 sites 


Track length 


— 


semi-infinite 


Track site spacing 




5.0nm x 5.0nm 


Initial set of product sites 


P 





Effective substrate binding rate 




1.0 x 10 3 s _1 


Effective product binding rate 


kp 


1.0 x 10 3 s _1 


Substrate dissociation rate 


k s 


0.0s- 1 


Product dissociation rate 


kp 


1.0s- 1 


Catalysis rate 




< 1.0s- 1 


Temperature 


T 


300 K 


Force in — x direction 


f 


<4.0pN 


Largest simulated time 


^max 


< 1.0 x 10 7 s 



it represents the sole kinetic difference between substrates and 
products. We have fixed k$ = kp so that there is no attach- 
ment bias between substrates and products. An unattached 
leg will just as rapidly bind to a feasible substrate as to a fea- 
sible product. Once bound, a leg-product complex unbinds 
at rate kp = 1, and a leg-substrate complex unbinds at rate 
+^cat- We assume that substrate unbinding is much less 
probable than substrate catalysis so we let &<7 = 0. Then if 
& cat = 1, there is no residence time bias between substrates 
and products — the expected duration of a leg-product bind- 
ing is the same as that for a leg-substrate binding. While 
substrates are still converted into products, the kinetics of 
the walker attachment and detachment are identical for both 
species. Hence, a walker with & cat = 1 is equivalent to a walker 
moving over an all-product surface. But an all-product surface 
provides no chemical free energy, and so an all-product walker 
system must move diffusively. The walkers with & cat = 1 still 
release chemical energy when they catalyze a substrate, but 
the kinetics prevent them from utilizing that energy, and so 
they represent the no-energy baseline motion of walkers. In 
contrast, when & cat < 1 there is a residence time bias, where 
leg-substrate bindings are longer in duration than leg-product 
bindings, as the leg must wait until the relatively slow catal- 
ysis step completes before it can unbind. The only part of 
the walker kinetics which takes into account the chemical free 
energy released in substrate catalysis is the assumption of irre- 
versibility in the enzymatic conversion from substrate to prod- 
uct. In enzyme kinetics there is some non-zero rate for the 
reverse of the catalytic process. However, if the Gibbs free 
energy (AG) drop from substrate to product is large enough, 
the reverse rate is so small it is for all practical purposes zero, 
and is omitted from the walker kinetics in our model. Thus, 
we vary the & cat parameter to control the residence time bias 
between visited and unvisited sites, and at & cat = 1 the motion 
of the walker is equivalent to the no-free-energy case. We do 
not directly incorporate AG into the model, as the kinetic val- 
ues of & cat and &p are more important to walker motion than 
AG, and any free energy change large enough to make the sub- 
strate modification effectively irreversible is sufficient to meet 
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FIG. 5: (color online) Simulation estimate of ( ||p(f) || 2 ) when / = 0. 
Walkers with fc cat = 1 move diffusively. Those with & cat < 1 move su- 
perdiffusively, but eventually use up their local supply of substrates 
and become ordinary diffusive. True transitions to diffusion will oc- 
cur above simulated time ^max = 

10 7 . 



the model assumptions. 



A. Walkers move superdiffusively 

Figure [5] shows the ensemble estimates (N = 1000) for 
MVRWs moving in the absence of a load force. Initially (be- 
low the characteristic timescale of l/k cat ) the walkers move 
subdiffusively. As expected, the k cat = 1 walkers never move 
faster than diffusion. However, as & cat is decreased, walkers 
initially move more slowly due to the the slower catalysis 
kinetics, but once sufficient time has passed, they move su- 
perdiffusively with a > 1. The smaller the value of k CSLt , the 
more superdiffusively the walkers move, with a approaching 
2 for the smallest k cat values. This superdiffusive behavior 
persists over several decades in time, during which the walk- 
ers are moving processively away from the origin in the direc- 
tion of unvisited sites. Because of this outward-directed bias, 
the walkers with k cat < 1 eventually overtake (in MSD) the 
kcat — 1 walkers given sufficient time. However, the ability 
to move superdiffusively depends on the local availability of 
the immobile substrate fuel, which is consumed as the walker 
moves over the track. Hence if a walker moves back over 
previously visited sites, it becomes starved for fuel. In these 
energy-devoid regions the walker can only move diffusively 
like the k cat = 1 walkers, and so superdiffusion must eventu- 
ally give way to regular diffusion, even for the smallest values 
of kcat. 

Figure [6] shows the number of sites catalyzed over time; its 
rate of change represents the average availability of substrate 
fuel. As long as the number of sites cleaved grows linearly 
with time, the walkers are receiving fuel at a constant rate and 
their motion is biased in the direction of new sites, which al- 
lows their constant fuel supply to be maintained. When N(t) 
becomes sub-linear the walkers begin to transition from su- 
perdiffusion to ordinary diffusion. 
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FIG. 6: (color online) Simulation estimate of (N(t)), the number of 
substrates catalyzed to products when / = 0. Since k$ = 0, this is 
equivalent to the number of distinct sites visited at time t. Walkers 
with & cat < 1 catalyze substrates at a nearly linear rate over many 
decades in time. This is necessary to maintain a constant supply of 
chemical energy to sustain superdiffusive motion. 



B. Walkers do work against a load 

To quantify the sensitivity of the walker's superdiffusive 
motion, we impose a constant conservative load force / on the 
walkers in the — x direction (Fig .14]). Figure [7] shows ensemble 
(N = 4000) estimates of ( ||p(f) ||^ under a range of forces for 
k cat = 1 and k cat = 0.01. Again, k cat = 1 (dashed lines) illus- 
trates the no-energy case and, as shown previously (Fig. [5]), 
walkers move diffusively without the influence of force. 

When / > 0, the random walk over products is biased in the 
— x direction. The constraints imposed by the surface prevent 
the walker from moving to the left of the origin, so the biased 
random walk will eventually reach an equilibrium position, 
after which the motion is stationary (a = 0). Indeed, this is 
seen for the k cat = 1 walkers, which never move faster than 
diffusion, and their MSD increases monotonically to the equi- 
librium value exactly as if they were undergoing constrained 
diffusion in a box (34). In contrast, when k caX < 1 we again see 
nearly ballistic motion for all walkers except those under the 
highest load forces / > 2.0 pN. Thus, even though the load 
force attempts to pull the walker body away from the sub- 
strate fuel, the long residence time for leg-substrate binding 
allows a few substrate-bound legs to resist the force and keep 
the walker in proximity to the substrate sites. Eventually as 
in the / = case, all walkers regardless of k cat will exhaust 
their local supply of substrates and will find themselves mov- 
ing over energy-devoid product sites, which ultimately brings 
them to the same equilibrium position as the & cat = 1 walkers 
(for a given /). 

The change in potential energy of the walkers as they move 
in opposition to the load force can be quantified by evaluating 
the ensemble estimate of the mean position of the walkers' 
bodies, (p(f))- We chose to set AE = when p x = 0, and then 
AE = fp x > for walkers to the right of the origin (Fig. [4]). 
Figure [8] shows the ensemble estimate of (AE(t)). As the load 
force is increased above 0, the walkers attain progressively 
higher potential energies, and their peak energies come earlier, 
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FIG. 7: (color online) Simulation estimate (N = 4000) of (||p(f)ll ) and 95% confidence bounds (shading) on a log-log scale. Reference lines 
are shown for ordinary diffusion (a = 1) and ballistic motion (a = 2). Walkers with & cat < 1 move superdiffusively, but when / > 0, they 
eventually slow down and return to the same equilibrium position as the & ca t = 1 walkers. 




FIG. 8: (color online) Simulation estimate (N = 4000) of (AE(t)) and 95% confidence bounds (shading) on a log-linear scale. Walkers with 
/ = always have AE = 0. Those with / > and & cat < 1 do significant amounts of work, reaching a peak energy before eventually coming 
to an equilibrium. This equilibrium depends only on / and not on k cat . 



as they need to move less distance to do the same amount of 
work. However, as the forces are increased beyond / = 2pN, 
the walkers are not able to move very far without being pulled 
backwards, away from their substrate fuel, and they achieve 
only modest values of AE. 



C. Peak work 

When / > all walkers eventually move to an equilibrium 
position with energy AE OQ (f). This value is greater than the 
initial energy, because the walkers begin out of equilibrium 
with only a single leg attached (Fig. [4]). The initial energy 
of the walker AEq(/) < because we measure p as the body's 
equilibrium position (B), which under any non-zero force will 
have p x < at the initial walker attachment location. How- 
ever, the kinetics of &p ^> &p lead to an equilibrium where 
legs are almost always attached to a site, and because all sites 
are to the right of the origin, the equilibrium position AE^ (/) 
will also necessarily be greater than AEo(f). Thus, to charac- 
terize the amount of useful work that a walker can do we take 
into account the equilibrium energy specific to each force. We 
define the peak work for force / as 

**(/)= max <AE(f;/))-A£o(/). (9) 

te[o,t max ] 

We estimate AEoo(f) as (AE(t max ;f)) for the k cat = 1 walker. 
Figure [9^ shows w* as force and & cat are varied. The k cat = 1 
walkers never have w* > 0, but the walkers with & cat < 1 can 
do significant work under moderate forces. We also show the 
values for p*(f) = max^ [0/max] (p x (t\ /)) - p~ '(/) in Fig. [5k 
Note that the walkers move significantly farther under small 
loads, although they do nearly the same work. 



V. DISCUSSION 

Multivalent random walkers are able to do work because 
they act as Brownian ratchets. The physical motion of the 
walker is the result of random thermally driven molecular mo- 
tions that are rectified by the constraints imposed by attached 
legs. Without any structural or conformational coupling, the 
independently operating legs are constrained only by their 
passive connection to a common body. The gaits with which 
MVRWs move are uncoordinated, unoriented, and acyclic, yet 
they can be designed to move nearly ballistically along tracks 
laid out in 2D. 
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FIG. 9: (color online) Simulation estimate of (a) peak work w*(f) 
and (b) peak x position /?£(/). Walkers with & cat G {1,0.01} are 
shown with an ensemble estimate of the mean using N = 4000 sam- 
ples; other £ cat are shown using N = 250 samples. The peak position 
for / = is shown as well, which is limited by the simulated time 
*max = 1.0 x 10 6 . In particular when / = 0, the & cat = 0.001 walk- 
ers are still moving superdiffusively at t = 1.0 x 10 6 , but are limited 
by their slower stepping kinetics. At longer times the & cat = 0.001 
walkers will achieve a peak position greater than those achieved by 
the larger & cat walkers. 



to form (Fig. [TO]). At the boundary between the product sea 
and unvisited substrates, the local substrate concentration gra- 
dient is in the -fx direction [44]. The emergence of spatial 
asymmetry in concentration makes it possible for an unori- 
ented, symmetric walker to develop a directional bias. At the 
boundary, a MVRW with & cat < 1 is biased in the + x direc- 
tion not because the legs are more likely to attach to substrates 
(&+ =kp), but because when they do attach to a substrate, they 
stay bound longer — there is an effective residence time bias. 

A walker with & cat < 1 is only directionally biased when 
near the boundary, in which case its legs irreversibly catalyze 
attached substrates to products, moving the boundary in the 
+ x direction as well. Thus, as long as a walker remains near 
the boundary, it is biased in the + x direction, and it moves the 
bias-inducing substrate concentration asymmetry along with 
it, which leads to persistent motion directed away from the 
origin. 

In ID without force, the motion at the boundary has been 
shown analytically to be ballistic fl4l . From results in 
Figs.|5]and[7]the motion in 2D (at the ensemble level) is nearly 
ballistic even when opposing small forces, implying that indi- 
vidual walkers near the boundary must also be moving nearly 
ballistically. 



A. Mechanism of superdiffusive motion 

The superdiffusive motion of walkers and its eventual decay 
to diffusion (/ = 0) or stationary equilibrium (f > 0) can be 
understood by noting that the only source of energy available 
to the walkers is present in the substrate molecules, which are 
a locally-limited, immobile resource. 

After the walker starts moving and catalyzing sites, a con- 
tiguous region of product sites we call the product sea begins 



B. The boundary and diffusive metastates 

The emergence of the boundary between the product sea 
and the unvisited substrates causes the walker to move su- 
perdiffusively, but eventually all walkers either move diffu- 
sively (/ = 0) or move to a stationary equilibrium distribution 
(f > 0). This behavior can be understood by decomposing 
the Markov process into two metastates: a boundary state (B) 



FIG. 10: (color online) The irreversible catalysis of substrates to 
products leads to the emergence of a spatial asymmetry in substrate 
concentration at the boundary between the contiguous product sea 
and the contiguous region of unvisited substrates. A walker with 
& ca t < 1 has a residence time bias where leg-substrate binding dura- 
tions are much longer than for leg-product bindings. Thus, over time, 
legs are more likely to be attached to local substrates than local prod- 
ucts, not because they seek out substrates, but because legs attached 
to products quickly detach. Hence, walkers are effectively driven in 
the direction of greatest local substrate density, and near the bound- 
ary this is always in the + x direction. The irreversibility of substrate 
catalysis means the boundary itself also moves in the + x direction, 
causing walkers near the boundary to move ballistically away from 
the origin. 



wherein the walker is attached to substrates near the bound- 
ary, and a diffusive state (D) wherein the walker moves over 
the energy-devoid product sea (Fig.fTT]). 

When the walker is in the B state it moves ballistically in 
the + x direction, but when it is in the D state it has no direc- 
tional orientation, and it moves by ordinary unbiased diffusion 
for / = 0, or by — x-biased diffusion when / > 0. Figure 12 
shows three typical traces of the position of individual walk- 
ers under no force, where B and D periods have been shaded 
to show the alternation between states and the distinction be- 
tween the ballistic and diffusive motion. 

The probability of a walker leaving the B state by mov- 
ing sufficiently far in the — x direction is independent of the 
absolute position of the boundary. Thus, the B metastate is 
Markovian since the transition rate to the D metastate is inde- 
pendent of how long the walker has been moving or the size 
of the product sea. As & cat is decreased, the duration of leg- 
substrate bindings relative to leg-product bindings increases 
and the walker is less likely to simultaneously detach from all 
boundary substrates and leave the B state. Thus, lower & cat 
values result in more persistent ballistic motion over longer 



durations, but at smaller velocities (Fig. 12). 

In contrast the D metastate is non-Markovian. The dura- 
tion of a D period depends on the size of the product sea, and 
hence this duration grows as the walker catalyzes sites. In the 
case where / = 0, the time is quadratically dependent on the 
size of the product sea, but when / > this dependence be- 
comes exponential, and for sufficient forces and sufficiently 
sized product seas, the probability of returning to the bound- 
ary once departed a significant distance becomes effectively 
0. Hence, the duration of ^-periods is constant in time, but 
the duration of £>-periods grows. Eventually walkers spend 
nearly all their time moving over products in the D state, and 
so approach the same equilibrium distribution as the k cat = 1 
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FIG. 11: (color online) The walker moves between boundary (B) 
and diffusive (D) metastates. The walker moves ballistically in the 
direction of local substrate gradient when in the B state, but moves 
diffusively over previously visited sites in the D state. The walker 
initially spends most of its time in the B state, consuming substrate 
fuel, however as the product sea grows, the time to exit the D state 
increases, leading to asymptotically diffusive motion in the absence 
of force and equilibrium stationary motion in the presence of force. 



walkers, as seen in Figs. [7] and [5] 



C. Comparison with kinesin 

Multivalent random walkers accomplish processive, su- 
perdiffusive transport through significantly different mecha- 
nisms from known natural molecular motors, such as kinesin. 

• Kinesins move over regular, oriented polymeric tracks, 
while MVRWs move over arbitrary, unoriented 2D sur- 
faces of substrate tracks. The MVRW tracks can be 
heterogeneous or homogeneous, structured or unstruc- 
tured. 

• Each kinesin head has an oriented binding site for 
the track and a separate binding site for the substrate 
fuel. These two binding sites are kinetically coupled 
by complex conformational interactions Q. Further- 
more, the actions of the two kinesin heads are chemo- 
mechanically coupled by long-range conformational 
changes (7), whereby the kinetics of a single head dis- 
sociating depends on the chemical state of the other 
head l35l . In contrast, the legs of a MVRW are un- 
oriented and uncoupled. The conformational state of an 
individual MVRW leg has no influence on the chemical 
or mechanical actions of the other legs. 
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FIG. 12: (color online) Typical traces of p x for a MVRW with / = 
for three & ca t values. The traces are shaded blue when the walker is in 
the B metastate, and red when it is in the D metastate. Walkers with 
smaller & cat have longer B periods, but smaller velocity. The duration 
of 5-periods is independent of time, but the duration of Z)-periods 
grows with the size of the product sea, and consequently increases 
over time. Thus, at short times the walker is more likely to be in the 
B state, but at longer times is more likely to be in the D state. 



• Kinesin always moves with a rigid hand-over-hand 
walking gait where the leading head alternates at each 
step. The chemomechanical coupling ensures that the 
trailing head understands its orientation with respect to 
the track and leading head, making it more likely to un- 
bind and move than the leading head 0. A MVRW, 
on the other hand, has no preferred gait and legs are 
unaware of their relative positions or orientations. The 
leg motion is independent and uncoordinated, restricted 
only by the finite length of legs, and their connection to 
a common body. 



D. Dissociation 

There is a non-zero probability for a walker to detach from 
the track if A: — 1 legs are simultaneously in the detached state, 
and the next action chosen is for the remaining leg to detach. 
A walker with k detached legs is free to diffuse in solution, 
and cannot be ascribed a well-defined position with a discrete 
state Markov process. Hence, dissociation poses mathemati- 
cal difficulties for analyzing a non-ergodic motive process and 
comparing it with other mathematical models of anomalous 
diffusion. Ergodic models of kinesin can simultaneously ana- 
lyze motion and dissociation, because the transport character- 
istics and dissociation probabilities can be understood inde- 
pendently by studying a single motor cycle 13611571 . MVRWs, 
being non-ergodic, have transport and dissociation probabili- 
ties that depend on the current state of the local chemical sites, 
and cannot be analyzed with similar techniques. 

One approach to dealing with dissociation in non-ergodic 
walker models is to have a single absorbing dissociated state 
to which all walkers will eventually go and never return. This 
state is then the single equilibrium state of the system, and 
analysis is done on the remaining walkers. However, analyz- 
ing MSD becomes challenging because at any t > there is 
necessarily some non-zero proportion of walkers in the dis- 
sociated state. Ensemble MSD is no longer well-defined, as 
we cannot ascribe a position to dissociated walkers. Instead 
of this approach, we implement a hopping rule, whereby a 
walker with k—l legs whose next KMC chosen transition is 
to detach its one remaining leg is prevented from diffusing 
away from its dissociation location. It is temporarily held in 
place until a leg attaches to a local feasible site. The net ef- 
fect is a hop from one site to another and is implemented as a 
single KMC step, so that the position of the walker is always 
well defined. 

For any finite k$ and kp rates, it is possible for walkers to 
temporarily dissociate via a hopping event. In practice, how- 
ever, when the walker has sufficiently many legs, the on-rates 
are sufficiently fast, the legs are long, and the substrates are 
densely spaced, the probability of dissociation is low. Over 
the course of the simulations shown in Figs. [7] and [8] only four 
out of 56000 walkers with / < 3.0pN and 100 out of 16000 
walkers with / > 3.0pN experienced any hopping event. 



E. Effect of variation of number of legs and leg length 



In each of these areas of difference between kinesin and 
multivalent random walkers, the mechanisms employed by ki- 
nesin are more complicated to engineer from chemical build- 
ing blocks, but lead to more efficient and processive motion. 
Yet, without any of these complex mechanisms MVRWs can 
still transform chemical energy into mechanical work, which 
implies that none of these properties are actually necessary 
for a molecular walker to act as a translational molecular mo- 
tor. Thus the value of the multivalent random walker model 
derives from its simplicity — it has been stripped of all extra- 
neous properties leaving only those core features which are 
necessary for non-diffusive motion and mechanical work. 



As summarized in Table [TJ our results focus on four-legged 
walkers with leg length £ = 12.5 nm, which is 2.5 times the 
5.0nm substrate spacing distance. Both the leg length and the 
number of legs can be freely varied. However, there are sensi- 
ble ranges for these parameters, outside of which the motion 
of the walkers is not as processive, or is exceedingly slow. To 
be efficient molecular transport devices, walkers must simul- 
taneously avoid dissociation, resist the effect of forces, and 
remain attached to substrates near the boundary. 

First, consider the number of legs, which is varied in the 
range 2 < k < 5 in Fig. 13 For the residence time bias to 

we require k > 2 [14]. With few 



in Fig. 
lead to a directional bias 




10 1 10 2 10 3 
Time t (s) 



FIG. 13: (color online) Simulation estimates (N = 400) showing 
the effect of the number of walker legs (k) on walker motion when 
/ = 0. The MSD is shown with 95% confidence intervals in shading. 
Walkers with more legs move with smaller diffusion constant when 



1.0 = 



= & P and there is no residence time bias. However, when 



&cat = 0.01, the walkers with more legs experience a stronger di- 
rectional bias towards the local substrate concentration gradient and 
hence move superdiffusively over longer times and distances. Of the 
configurations studied, walkers with k = 5 legs and k cat = 0.01 even- 
tually achieve the greatest mean squared displacement. The black 
lines show the case k = 4, corresponding to walkers in Fig. [5] 



legs (k = 2), walkers are more likely to have all of their legs 
detached simultaneously and undergo a hopping step. As the 
number of legs is increased this probability drops exponen- 
tially, as each leg's probability of detachment is approximately 
independent. Walkers with more legs also tend to move more 
superdiffusively and processively when in the B state, as they 
have a higher probability that at least one leg remains attached 
to a substrate at the boundary. However, walkers with many 
legs have a significantly smaller diffusion constant. Hence, 
k = 4 was chosen as a reasonable compromise value that pre- 
vents dissociation, maintains a strong tendency to remain on 
the boundary, and moves appreciably fast. 

The leg length £ must be considered in relation to the sub- 
strate spacing. Together these parameters determine the av- 
erage number of feasible sites available for attachment. Ide- 
ally, the substrate spacing would be made as small as possi- 
ble, but it is constrained by limits on the sizes of molecules 
and how closely substrates can be arrayed on a surface. We 
chose 5.0nm as a reasonable lower limit on this spacing, as 
it approximates the density of DNA substrates arrayed on a 
DNA origami [38 ] surface, such as the substrate density used 
in molecular spider experiments ifTTTl . 

Figure [14] shows the effect of varying the leg length for 4- 
legged walkers while keeping the substrate spacing constant. 
We find that if legs are too short (£ < 5.0 nm), the number of 
feasible sites is too small to maintain a superdiffusive effect. 
For leg lengths i > 7.5 nm, which is 1.5 times the substrate 
spacing of 5.0 nm, there is little qualitative difference in the 
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FIG. 14: (color online) Simulation results (N = 400) showing the ef- 
fect (at / = 0) of varying the leg length 5.0 nm <£< 15.0 nm, while 
the number of legs is fixed at k = 4 and the substrate spacing is fixed 
at 5.0 nm x 5.0 nm. The effect of £ on MSD is shown with shading 
indicating the 95% confidence interval for the mean. The effect of 
changing the leg length is essentially manifested as a change in the 
diffusion constant, but not in the qualitative characteristics of the su- 
perdiffusive motion for the k cat = 0.01 walkers. The exception is for 
the very short leg length £ = 5.0 nm, where the average number of 
feasible sites becomes so small that walkers lose their superdiffusive 
transport behavior 



walker motion, although longer legs do lead to a faster diffu- 
sion constant in the absence of force. Under load, however, 
leg length and substrate spacing should both be minimized to 
maximize the peak work and displacement of walkers. Longer 
legs allow a larger feasible region T, leading to a larger bias in 
B under any non-zero load. This in turn makes it more likely 
for long-legged walkers to move backwards. We found that 
a leg length of approximately 2.5 times the substrate spacing 
provides a good balance between dissociation and processiv- 
ity, although a full analysis of this relationship is reserved for 
future study. 



F. Sensitivity to kinetic parameters 

The MVRW model has too many independently variable ki- 
netic parameters to simultaneously examine the effect of each 
of them on walker motion characteristics. We have chosen 
representative kinetics summarized in Table [I] to act as a ref- 
erence point. Clearly, any chemical implementation of the 
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We show the results of varying k s in Fig. 16 In other re 
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FIG. 15: (color online) Simulation results (N = 1000) showing the 
effect (at / = 0) of varying k£ on the MSD of walkers while kp 
is fixed at 1.0 x 10 3 s _1 . The black lines represent the case where 



kp 



^5 = l.Ox 10 3 , which is used in all other simulations results. 



When &p = k£, there is no attachment preference for a substrate over 
a product. An unattached leg will just as rapidly bind to a feasible 
substrate as to a feasible product. However, as we take k$ < k£ , the 
walkers have an attachment bias to products, which should be ex- 
pected to reduce the time spent in the boundary (B) metastate, and 
therefore lead to less pronounced superdiffusive behavior. These re- 
sults show that the MSD is robust to moderate changes in the on- 
rates, and even for k£ = k£/ 10, there is an appreciable superdiffusive 
effect when & cat = 0.01. However taking k$ = k^ /100 overwhelms 
the residence time bias of the walkers in the B state and prevents any 
superdiffusive motion. 



multivalent random walker model (e.g., molecular spiders) 
will have potentially very different (and likely much faster) 
rates than those we have chosen. However, it is not our pur- 
pose to model a specific chemical implementation. Instead, 
we show that the qualitative characteristics of superdiffusive 
walker motion persist over a wide range of kinetic values, as 
long as the residence time bias between visited and unvisited 
sites leads to an effective bias in the direction of local substrate 
gradient. 

In Fig. [15] we show that the superdiffusive behavior as quan- 
tified by MSD persists over an order of magnitude in variation 
of the k$ rate. Indeed, even if the walker is biased 10:1 in at- 
tachment preference to products over substrates, the residence 
time bias of 100:1 of substrate to product binding duration is 
still sufficient to achieve a superdiffusive scaling of MSD over 
several decades in time. This robustness even to large changes 
in attachment rates allows us to be confident that superdif- 
fusive behavior is a pervasive feature of multivalent random 
walker systems and is not critically dependent on our particu- 
lar choice of attachment rates. 



suits we have assumed that k s = 0, which is reasonable as 
this rate is likely to be much slower than kp or & cat for any 
practical enzymatic implementation of a multivalent random 
walker. Figure [16] shows that indeed the superdiffusive be- 
havior is robust to changes in as long as it remains sig- 
nificantly slower than kp and & cat . However, setting k$ = kp 
eliminates any superdiffusive effect, as there is no longer a 
residence time bias between substrates and products, and the 
motion of the walker near the boundary is no longer biased in 
the direction of the local substrate concentration gradient. 



G. Effect of forces on dissociation reactions 

In describing our model we show how forces affect the bi- 
molecular association rates through Eq. [6] However, apply- 
ing a load force to walkers should also affect the kinetics of 
the unimolecular dissociation events. From a high-level view, 
the kinetics of unimolecular reactions depend on a molecule 
having enough internal energy to surmount some reaction en- 
ergy barrier Uq. Thus, rate laws follow the Arrhenius formula, 
k(T) oc exp (—Uo/kBT). The effect of a force / applied to the 
molecule is a mean change in energy of AC//, and the rate is 
modified to 



k(T) = VGxp((AU f -Uo)/k B T). 



(10) 



The value of the constant v and the relationship of AU f with 
force / depend on the specific internal chemistry of the leg 
tethers, enzymes, and substrates [ 39 -41 ], the details of which 
are beyond the scope of our coarse-grained walker model. We 
surmise that the effect of small forces is a slight increase in 
fc7 and &p , although this change would not be uniform over 
all legs, as those attached to sites further in the -fx direc- 
tion will oppose more of the load force on average than other 
small increases in k$ do not qual- 



16 



sites. Based on Fig. 
itatively change the motive properties of the walker with re- 
gard to MSD, except when the forces are large enough so that 
+^cat > &p , which eliminates the residence-time bias and 
all superdiffusive motion. Additionally, small increases in kp 
actually lead to an increase in the MSD, as they increase the 
residency time bias. Quantitative analysis of variation in kp 
is beyond the scope of this paper, but is available along with 
other supplementary results for the MVRW model l42l . Over- 
all, these results show that even though the present formula- 
tion of the MVRW model does not describe the effect of force 
on dissociation rates, we expect an extension of the model in- 
cluding these rates to predict similar superdiffusive behaviors, 
as long as the forces and corresponding rate changes are small. 
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FIG. 16: (color online) Simulation estimates (N = 1000) of the mean squared displacement of the walkers as k s is varied. Shading shows 95% 
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shows the case k$ = 0. The value of k$ determines the rate of detachment without enzymatic conversion of the site to a product. As long as 
+ ^cat < &p » there remains a residence time bias, and our results show that walkers with & cat = 0.01 move superdiffusivly for k$ < kp = 1. 
This shows that the qualitative behavior of the walkers is unchanged for small variations in k$ , and the choice of k$ = in the model is 
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eliminated, as there is no longer an effective residence time bias between visited and unvisited sites. 
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